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erbium in YAR By the investigation of the 
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two crystalline directions, the Er was 
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the impurity atoms used for the applica- 
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Introduction 



GRID (Gamma Ray Induced Doppler broadening) 
spectroscopy is based on the analysis of a Doppler 
broadened secondary 7 line which is emitted from 
a nucleus highly excited due to neutron capture and 
recoiled after emission of a primary 7 quantum. The 
Doppler broadened y line contains information about 
the lifetime of the excited nuclear states as well as 
the stopping cross section of the recoiling nucleus in 
the solid [1.2]. The method has been mainly used in 
polycrystalline samples to determine the lifetime of 
nuclear levels and to investigate the interaction of low 
energy ions with atoms of solids [3,4]. 

If the recoiling nucleus is embedded in a single 
crystal its slowing down is non-isotropic. It depends 
on the surrounding of the recoiling atom and on the 
crystalline direction. Molecular dynamics (MD) simula- 
tions and experiments on single crystals have shown that 



the energy distribution of the Doppler broadened 7 line 
depends on the orientation of the crystal in relation to 
the direction of observation [5-8]. Consequently, GRID 
experiments in combination with MD simulations can 
give information about lattice positions of foreign atoms 
as well as about the structural disorder surrounding the 
recoiling atom. 

In the present paper some conditions necessary for 
the application of GRID spectroscopy are discussed. As 
an example, a first test of GRID spectroscopy for lattice 
site localization of Er in yttrium-aluminium-perovskite 
(YAP, YAIO3) is shown. 

2. Experimental 

The applicability of the GRID spectroscopy for 
materials investigation is restricted to elements having 
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an adequate 7 cascade with a sufficiently high intensity 
after neutron capture, and the lifetime of the inter- 
mediate energy level emitting the secondary 7 radiation 
should be in such an order of magnitude that the recoil- 
ing nucleus may reach the next, or the next but one 
neighbouring atom. The Doppler broadening of the 
secondary 7 line of the recoiling nucleus is typically in 
the order of t^E/E ~ 10 * to 10"^ depending on the 
primary 7 emission and the nuclear mass. And finally, 
the probability of the 7 emission is determined by the 
thermal neutron flux for the initial capture process and 
the number of recoiling nuclei per target. The resulting 
intensity has to be compared with the resolution and the 
luminosity of the 7 spectrometer to get an estimate 
about the possibility of detecting any structures of the 
Doppler broadened line profile. 

The first GRID experiments on foreign atoms in 
single crystals were carried out at the Institute 
Laue-Langevin in Grenoble. The crystals were mounted 
near the reactor core where a flux of thermal neutrons 
of 10''' cm"^ s ' is maintained. The 7 line profile was 
monitored along two crystalline directions with the 
double flat crystal spectrometer GAMS4 [9] having an 
energy resolution in the order of bJL/E ~ 10 '' to 10 * 
depending on the detected 7 line. 



3. Results and Discussion 

3.2 Experimental Results on Erbium in YAP 

For the first test of the GRID spectroscopy for foreign 
atom localization we chose Erbium which is an interest- 
ing element for both nuclear and solid state physics. 
We analyzed the 1649 keV line of "'Er (1913 keV -^ 
264 keV) along the [001] and [201] directions. Pretests 
showed that the 1649 keV transition exhibits a Doppler 
broadening and that the 1913 keV level has a half-life- 
time shorter than 70 fs. The parent "'Er has naturally an 
isotopic abundance percentage of /iso = 20 % and a 
neutron capture cross section of cr = 700 b. However, 
due to the low probability of the chosen 7 transition 
(5X10"^ per neutron capture), we had to take a host 
lattice where a high Erbium concentration can be 
incorporated. Therefore, we took the laser material 
yttrium aluminium perovskite (YAIO3, YAP with a 
concentration of natural Erbium in the order of « = 
10^^ cm"^ The size of the targets was 20 mmX 
18 mmX2 mm. 

In Fig. 1, the experimental results for the [201] direc- 
tion along the Y-Al row (corresponding to the [111] 
direction in pseudo-cubic description) are compared 
with the calculated line profiles assuming that Erbium 
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Fig. 1. Measured energy distribution of the 1649 keV 7 line of "' Er in YAP along 
the [201] direction, the solid and dotted lines represent the calculated line profile 
for Er on the Y and Al sites, respectively. For comparison with the measured line, 
the simulated profiles of Fig. 2 and the response function (dashed line of the 
7 spectrometer are superimposed. 
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occupies either the yttrium site or the aluminium site. 
The calculations take into account the energy resolution 
of the spectrometer by superimposing the response 
function and the simulated line profiles of Fig. 2. Most 
of the structures of the line profiles in Fig. 2 are covered 
by the energy resolution which is low compared to the 
details of the line profiles, despite of the absence of 
pronounced structures in the measured line profile, we 
should show in Ref. [10] the the MD-simulations fit the 
experiment for both measured orientations if the erbium 
is assumed to occupy the Y site and the half-lifetime Tm 
of the 1913 keV level is about 20 fs. 

From the GRID measurements on Erbium in YAP we 
deduce the minimum 7 intensity required for GAMS4 
and the future GAMS5 operating with two curved crys- 
tals in DuMond geometry. The minimum intensities 
listed in Table 1 are estimated under the assumption that 
the measurement time should not be reasonably longer 
and the statistical error should be similar. The values 
obtained for GAMS5 are estimated with the lO'' times 
higher solid angle and the different target sizes [11]. 

For the application of GRID on impurity atoms, the 
minimum required impurity density in the crystal is of 



practical interest. The measured 7 intensity is propor- 
tional to the 7 intensity emitted from the target volume. 



oc o-XL^X«X/„„ = cr„tfXnXL, 



(1) 



o" represents the neutron capture cross section ly^ the 
intensity of the secondary 7 emission, n the density of 
the incorporated element, and /iso the isotopic abundance 
of the reacting isotope. The product of the neutron cap- 
ture cross section crand the intensity ly^ of the secondary 
7 transition is the effective cross section a^tf. 

The 7 spectrometer GAMS5 is sensitive to only a thin 
layer with a thickness of 20 jjim. Therefore, it can be 
applied to layers produced by MBE (Molecular Beam 
Epitaxy) or to impurity atoms incorporated by ion im- 
plantation, for these thin layers, the layer density D is a 
more common quantity, and the measured 7 intensity is 
proportional to 



a (tXL^XDXL, 



(2) 



Taking into account that the maximum density of 
an element in crystals is about n = 5Xl0^^cm"^ the 
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Fig. 2. Simulated 1649 keV line of Er in YAP along the [201] direction for Er 
occupying an Al site or a Y site, respectively. 



Table 1. Minimum intensities of the 7 cascades for the high resolution 7 spectrometer GAMS4 at ILL and the expected limit for GAMS5 in 
DuMond geometry 
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effective cross section of the y emission should be 
higher than 0.02 b when GAMS4 is used and can be 
reduced to 5X10"^ b for GAMS5 in future. In case of 
nuclei having transitions with higher effective cross 
sections a^ti, its concentration n in the crystal could 
be reduced until the minimum intensity for the used 
spectrometer is reached. 

3.2 Selection of Impurity Atoms 

Nuclei suitable for the GRID spectroscopy should 
show 7 cascades with a minimum intensity as given 
above. Additionally, the lifetime of the intermediate 
level should be in such an order that the recoiling 
nucleus interacts with the surrounding crystal and influ- 
ences the structure of the Doppler broadened 7 line. 

For getting a more precise idea about the influence of 
the lifetime on the Doppler line-profile, we carried out 
various MD simulations on different crystal lattices and 
elements. The results of these simulations shows that for 
the majority of lattice impurity combinations, the influ- 
ence of the impurity location on the Doppler broadening 
is most pronounced, if the mean range of the recoiling 
nucleus is in the order of 2 A to 3 A. This mean range 
corresponds to the atomic bonding length and conse- 
quently to the distance to the next neighboring atoms 
where the recoiling nucleus is scattered. 

For mean ranges shorter than 2 A, the majority of the 
recoiling nuclei decay before interacting with the next 
neighboring atoms and the lattice symmetry of the host 
crystal does not affect the Doppler line. 



For mean ranges larger than 3 A, the stopping of the 
recoiling nucleus by the next but one neighboring atoms 
becomes more and more important. The interaction with 
these atoms causes a further slowing down and reduces 
the width of the Doppler broadened line. Characteristic 
structures due to the initial location of the recoiling 
atom are more and more wiped out [10]. 

However, for a pre-selection of suitable nuclei, we 
looked for a simple approximation which takes into 
account only the nuclear parameters of the recoiling 
nucleus. We found that the product of the initial velocity 
of the recoiling atom and the half-lifetime of the inter- 
mediate level can be used, subsequently called effective 
mean range Sdf, 



■5eff — VoXTi/2, 



with 



Vo = Ey'/mc . 



(3) 



The initial velocity Vo is given by the energy E^' of 
the primary 7 transition, the mass m of the recoiling 
nucleus, and the light velocity [3]. 

Compared to the mean ranges obtained from the MD 
simulations s^ff is independent of the used elements of 
the crystal and it can be calculated easily. For instance. 
Fig. 3 shows the effective mean range of Erbium as a 
function of the half-lifetime of the 1913 keV level. 
Additionally, the mean ranges from simulations of 
erbium in three different crystals are drawn. The 
comparison of the ranges shows, that for the most 
interesting range up to 3 A, s^tf is sufficiently close to 
the simulated mean range to be used as a criterion for 
the pre-selection. 
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Fig. 3. Mean range of the recoiling erbium as a function of the half lifetime of the 
1913 keV level. The simulated mean ranges are depending on the crystal where the 
erbium is incorporated. The straight line indicates the effective mean range s^ff. 
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For ^cff larger than 3 A, the difference between the 
mean range from MD simulation and s^t! increases 
depending on the used crystal and the nucleus. A simple 
estimation is no longer possible. 

Taking into account the restrictions concerning the 
intensity of the 7 cascade [Table 1, Eqs. (1), (2)] and the 
lifetimes of the intermediate level [(Eq. (3)], ail known 
7 transitions listed in the database ENSDF published by 
the Lawrence Berkeley National Laboratory were 
analyzed. Additionally, the own results on Erbium 
"''Er(n,7) from the GRID experiment are included [10]. 

In Fig. 4 the neutron capture cross sections and the 
intensities of the secondary 7 transitions belonging to 



the suitable 7 cascades of the nuclei are indicated. Some 
nuclei like Ti show several useful 7 cascades. For a 
better overview, only the energy of the secondary 
7 transition belong to the most promising 7 cascades are 
indicated in the figure. 

Figure 5 gives an overview about the effective ranges 
ieff of the nuclei before emitting the 7 lines of Fig. 4. Of 
all known decays, only Si, Ti, Cr, Ba, and Er show 
7 cascades leading to effective ranges close to the range 
from 2 A to 3 A. All other nuclei like CI or Fe show 
longer effective ranges with the consequence to be 
examined more carefully before application in GRID 
experiments. The effective range of "''*Er fits to that 
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Fig. 4. Neutron capture cross sections and effective intensities of the strongest 
secondary y transitions belonging to the y cascades which lead to effective ranges 
in the range of 1.5 A to 10 A. The energy of the most promising secondary 
7 emissions are indicated in keV. 



Fig. 5. Effective ranges of neuclei before emitting the secondary y lines shown in Fig. 4. 
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required for the GRID experiment and, despite of its 
low probability, the effective intensity of the 1649 keV 
transition is one of the strongest which could be found. 
Up to now, no suitable y transition stronger than 
o"eff = 4 b is found in the nuclear data sheets. This re- 
stricts the application of GRID to foreign atoms in solids 
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GAMS5. Especially for heavy nuclei, many y cascades 
and lifetimes are not sufficiently known yet. This gives 
the chance, that the list of suitable nuclei for application 
of GRID on impurity atoms in crystals will be extended. 

4. Summary 

MD simulations and first experiments on Er in YAP 
confirm that GRID experiments can be used to localize 
foreign atoms in crystals. Due to the nuclear parameters 
and the experimental conditions, the method can be 
applied only to a very low number of impurity atoms. 
With the present high resolution 7 spectrometers 
GAMS4 at ILL, the concentration of the impurities 
should be higher than 2.5 X 10'" cml With the GAMS5 
spectrometer, which is promised to be available in the 
near future, a minimum density of 6X 10" cm"^ will be 
possible. 
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